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THE FLEXIBLE EECTANGOLAR WING IN ROLL AT SUEERSONIC FLIGHT STEEDS 
By Warren A. Tucker and Robert L. Nelson 


SUMMARY 


The problem of the loss In rolling effectiveness at supersonic flight 
speeds is considered for the case of the rectangular wing. Equations are 
obtained from which can / be calculated either the loss in rate of roll due 
to flexibility op the torsional stiffness required to maintain a given 
rate of roll, for two assumed variations of the spanwise distribution of 
torsional stiffness. 

A computation form and figures are provided so that calculations can 
be made without reference to the details of the analysis. 


INTRODUCTION 


The problem of the flexible rolling wing in subsonic flow haB been 
dwelt upon by many investigators, who have devised numerous methods of 
solution. References 1 and 2 were found valuable in the preparation of 
the present paper; reference 1 notes several other useful papers. So far 
as the domain of supersonic speeds is concerned, however, little work 
has been done on the problem. 

The fundamental relations between the aerodynamic and the elastic 
forces are, of course, valid far any flight Bpeed. The aerodynamic forces 
themselves, however, are critically dependent, on the flight speed because 
of the vast difference in flow behavior between the subsonic and the 
supersonic regions. Any solution of the aeroelastic problem at supersonic 
speeds, then, must consider aerodynamic forces calculated by the various 
methods peculiar to the supersonic zone. 

As a first step toward a solution of the problem, a treatment is given 
of the relatively simple case of loss in rolling effectiveness of the 
rectangular wing with constant— chord ailerons extending inboard from the 
wing tips. (See fig. 1.) The methods used are in general based on the 
linearized equations for supersonic flow. Due consideration is given to 
three-dimensional effects. The torsional stiffness of the wing is assumed 
to vary either linearly with the distance outboard from the wing center 
line or as the cube of this distance. 
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Because of the relative complexity of the equations involved in the 
analysis^ a computation form and figures are provided from which calculations 
can he made -without reference to the details of the analysis. 


SYMBOLS 



Po 

1 



5 


e 

6 r 

pb/2Y 


aspect ratio (b/c or b 2 /s) 

wing span 

total aileron span 

wing chord 

aileron chord 

wing area (be) 

semispan (b/2) 

distance out from wing center line (see fig. l) 
free— stream Mach number 


ratio of specific heat at constant pressure to specific 
heat at constant volume (1.40 for air) 

free— stream static pressure 

free-etream dynamic pressure 

lifting pressure 
lifting-pressure coefficient (P/q.) 

additional, two-dimensional lifting-pressure coefficient 
due to aileron deflection (^8/0) 

deflection of one aileron 

angle of twist at any spanwise station 

angle of twist at reference station (midspan of aileron) 

wing-tip helix angle 
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p rolling angular velocity 

Y free-etream velocity 


l 

°l 



T 


m 

GJ 


rolling moment 

rolling-moment coefficient (l/qSb) 
rol lln g-mament-effectivenesB coefficient for two ailer- 

“ (I 1 ) 

damping-in-roll coefficient 
rolling-moment— loss coefficient 




torque at any spanwise station 



pitching moment about midchard at any BpanwiBe station 


torsional— stiffness parameter 



m 0 torsional— stiffness parameter (T/0) 

m0 r torsional-stiffness parameter at reference station 

(midspan of aileron) 

k • constant of proportionality for stiffness distribution 

v ±) v 2 . constants of proportionality for twist distribution 

F tr aili n g-edge angle reduction factor 

f disturbance— velocity potential 

Cartesian coordi na tes parallel and normal, respectively, 

to free-etream direction (see fig. 13(a)), for field 
points 

£,T) Cartesian coordinates parallel and normal, respectively 

to free-etream direction (see fig. 13(a)), for source 
points 



4 


NACA TIT No . 1769 


^¥-?(¥) 2 

Subscripts: 

B rigid 

0 sea level 

1,2,3 used with 0 and Cp to indicate various regions 

of wing 

All angles are in radians, unless otherwise specified. 

ANALYSIS 

Tn a steady roll, the resultant rolling moment act ing on the wing 
is zero. This statement can he written in equation form as 

- (-S%) If - 0%/r = 0 (1) 


where PC-^ is the roll ing-moment— effectiveness coefficient of the ailerons, 
is the damping-in-roll coefficient, and is the rolling-moment— 

Iobb coefficient resulting from the flexibility of the wing. The quantity P 
is ret ain ed merely for convenience in the subsequent analysis. The positive 
sign is given to the quantities PC^, — pCx^, 31111 0 C Z 0 • 

The following equations can be written from equation (l) : 


pb/2V 

5 






( 2 ) 
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p-b/27 i PC ^ r 9 r 
(pb/2V) E 3C2 5 8 


Equations (l) to (3) are applicable to flight at either subsonic or 
supersonic speeds. The problem, then, is simply to determine the values 
of the various quantities at supersonic speeds. For a given wing-aileron 
configuration, equations for pC^g anji H can 130 found from refei>- 


ences 3 and h, respectively. The remainder of the analysis is devoted to 
the derivation of expressions for the rolling-mament— loss coef- 
ficient and the twist parameter 0 r /8. Those readers who are 

more interested in the application than in the derivation are advised to 
turn to the section entitled "Computation Procedure and Discussion." 


Torque A knowledge of the spanwise distribution of twist 9 1b 

necessary for the determination of either pC-i or 0 r /5. The twist 

e r 

at any section is found from considerations of the aerodynamic torque and 
the elastic forces acting on the wing. The aerodynamic torque is assumed 
to result solely from the moment about the wing midchord line of the 
pressure resulting from aileron deflection. This assumption implies that 
the centers of pressure of the lift forces resulting from rolling and from 
twist are so near the elastic axis as to contribute only a negligibly mnnn 
twisting moment. Because the elastic axiB of a rectangular wing designed 
for supersonic speeds would be expected to lie near the midchord line 
(much as in the case of a subsonic wing, where the elastic axis wnfl the 
line of aerodynamic centers are not too far separated) anrl because the 
centers of pressure of the farces resulting from rolling and from twist 
are near the midchard line, the assumption appears reasonable. In 
addition, the loads caused by aileron deflection are assumed to be 
independent of wing twist, and the aileron itBelf is assumed not to 
twist. An important result of .mak in g these assumptions is the elimination 
of any necessity for using a successive approximation method. 

The aerodynamic torque T at any spanwise staition y is simply the 
total moment acting outboard of the station. That is. 


T = 



(*) 


where m is the pitching moment about the midchord at any station gnfl 
is assumed to result solely from the pressures caused by aileron 
deflection. These pressures were taken from the equations given in 
figure 2(a) of reference 3. The resulting equations for the torque were 
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very cumbersome; therefore, the exact expression vas approx imat ed by two 
simple expressions, one representing a constant torque over the part of 
the span inboard of the aileron (the value of this torque was taken equal 
to the greatest value given by the exact expression) and the other repre- 
senting a torque decreasing linearly along the aileron span to zero at 
the wing tip. The approximate equation for the torque per unit aileron 
deflection inboard of. the aileron is 


T = B 
Cp ro qc3 p 


(5a) 


where 0 = ^ ^ 1 — and the approximate equation for the torque over 
the aileron span is 


(5h) 


where 1 — ~ 5; A comparison of the torque distribution obtained 

from the exact and approximate equations is shown in figure 2 for a 
particular configuration and for two values of P. The same degree of 
approximation can be expected for other configurations. The principal 
factor governing the closeness of approximation is the value of P; 
as P increases, the approximation approaches more closely the exact 
variation. 

Twist .— The twist 0 at any s panwise station y can be expressed 


T _ B _b_ A, _ 7 \ 
qc3 P *a V- a/ 

ir oo 


0 = 



(6) 


therefore, if an espreBsion can be found for d0/dy, the twist can be 
calculated. 

How the torBional stiffness of the wing at any spanwise station can 
be represented either by the quantity GJ or by the quantity mg, where 
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GJ 


_ T 


de/dy 


— T 
% = - 


If the expression for GJ is integrated^ then 


6 

T 


JO 


dy 

GJ 


or 


mg 


Jo 


dy 

GJ 


bo that If j for example ^ 


then 




— = 
mg 


and if 


then 


■^ = hy2 

GJ ^ 


-I- = £ y3 

mg 3 J 




These results are of use in the subsequent analysis. 
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Consider the parameter GJ. From. equation ( 7 ) 


de _ _T_ 
dy GJ 


( 11 ) 


Substitution of this expression in equation (6) makes possible the 
calculation of 9 at any point along the span. 

The twist was calculated for two assumed spanwise variations of 
torsional stiffness. The first i = k (or i is known to he 

representative of many small-scale models used in free — flight supersonic 
research. -The second : = ky 2 ^ar ^ « 7^ is typical of many subsonic 

fighter aircraft. (The degree of applicability to future supersonic 
aircraft can only be conjectured.) The procedure outlined can be used 
to obtain equations for other stiffness distributions. The following 

equations are the results of the twist calculations: For — =» k. 


kCp qbc- 


JLZ 

20 s 


(12a) 


where 0 — % 1 — ~ and 

s b - 


B 


kCp^qbc- 






(12b) 


where 1 — ^ ^ = 1. For ^ = ky 2 s 


kCp qb3 c 3 2 ¥ 


= 2?B (s. 


(13a) 


y t}_ 

where ^ 1 — — and 
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0 _ B 

l&pjblcl 9 6p ^ ' 

where 1 — ^ < — < 1. 

d - s ~ 

The twist may "be expressed more conveniently as the ratio of the 
twist 0 at any station to the twist 0 r at a reference station, taken 

at the midspan of the aileron ^ = 1 - 1 From equations (12) 

for X= k , 



_0_ = ^a/b y 

0 r Bp s 


where 0 ^ ^ ^ 1 — enfl 
- s ~ h 



where 1 — 


<£<1 

= B = X. 


From equations (13) far 



0 _ ^b a /b / y \3 

0 r W 


where 0 ^ — = 1 
s 


h a 

flUll 

h 


(l^a) 


(W 


( 15 a) 
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9_ = 1_ 

9-p Eg 


k ?) 3 - at - i 1 



(15b) 


where 


_a < y < 
h s 


Although these expressions for the twist are not complicated, the 
use of a single twist equation which applies over the entire span results 
in a considerable saving of effort when evaluating the rolling-moment 
loss due to twist. For this reason, equations (llf) and (15) were 

approximated by the following expressions: For ■— = k (or — y ), 

CrJ \ JOQ / 


e 

e T 




8 - 

b 


(16a) 


where 


oSiS 

8 


1 


and for 




9 

9 r 



(16b) 


Y 

where 0 = = 1. These approximations are compared with eq ua tions (14) 

and (15) in figure 3 for two values of b a /b. 

Eolling-mament— Iobb coefficient 3Cb Q With the twist evaluated, 
^r 

there remains the problem of determining the rolling moment resulting 
from the twist. Equations (l6a) and (l6b) can be rewritten, respectively, 
as 
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from which the twist is seen to "be represented by two components, one 
linear with respect to y and the other parabolic with respect to y. 
The rolling-moment coefficient resulting from a linear twist has been 
evaluated in reference 4 and may be written as 



W 1 b 2 31.1 1.1 1 

V 2 3 2 W- 2 (PA) 2 S (PA)3 


(18) 


where w^/V is the coefficient of the y— term in the twist equation. The 

rolling-moment coefficient due to a parabolic twist is derived in the 
appendix and is 




w 2 
_2 b_ 1 

Y 4 2 





1 

(PA) 3 



(19) 


o 

where Wg/v is the coefficient of the y —term in the twist equation. 

If the expressions for w-^/V and W2 /Y from equations (17a) and (17b) 
are substituted in equations (l8) and (19) and if the derivative is 
taken with respect to 0 r , then after some simplification the following 

equations result for qj = ^ (° r — g 7 ) > 
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pc- 


0 r 8 - 5^ 


1 _ _1 + 1 * r 1 1 , 1 1 t 
2 pA 6 ( pA ) 2 2 



2 

3 


_1 

PA 



+ 



(20) 


and for 


_1_ 

GJ 



PC 




5 1 _ 1 1 _ 1 1 

3 (PA ) 2 2 (PA ) 3 8 (PA)^ 


(21) 


Twist parameter 8 r /5 .— The twist parameter may he evaluated directly 

in terms of nu , the torsional stiffness at the reference station 
°r 

(aileron midspan), from equations (12) and ( 13 ). The value of k is 
found from equations (9) and (10), and if the value of Cp is taken 
conservatively as k5/p (the linearized— theory result), then the following 

expressions are obtained: For — « y (or pir =^k 

Ulo \ LnJ 


9 -3 r„ 2BR 

r _ hf_ Pq. 1 

6 " ^ (W 3 K 1 - It) 


( 22 ) 
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and for 


1 


OC 

“e 


y 3 




ft 

6 


1)3 Pq 


BE 

2 


(PA) 3 

b 



( 23 ) 


COMPUTATION PROCEDURE AND DISCUSSION 


The solution of equations .( 2 ) and (3)' is most easily carried out 
with the aid of a computation form such as that shown in table I, which 
is largely self-explanatory. The -various quantities to he entered in 
the computation form are read from figures 4 to 11; in figures 6 to 11 
is noted the source of the data presented. The quantity PA, which is 
the abscissa for most of the figures, ie found from figure 5 for any 
value of Mach number and aspect ratio. The quantity Pq can be found 
either directly from table H for certain selected altitudes or from 
table II and figure 4 for intermediate altitudes. 

The determination of the rolling-naoment— effectiveness coefficient PC^ 

given in 

figure 6 represent linearized— theory values. The quantity F (see fig. 7 ) 
is a reduction factor which takes into account the finite trailing-edge 
angle of the aileron. The value of F was obtained by taking the value 
for the first two terms of the thickness reduction factor given in figure 1 
of reference 3. Some recent experimental evidence indicates that, 
presumably because of boundary— layer effects, the use of the first two 
terms gives results which are in closer agreement with experimental data 
than are results obtained by using the complete three terms. Although 
the factor F is strictly-' applicable only to flat— sided ailerons, it can 
be applied to other ailerons with little loss in over-al l accuracy by 
estimating an effective trailing— edge angle. 

An explicit solution for the reversal speed is difficult because of 
the manne r in which PA .enters the quantities involved. The reversal 
speed may be fo und by a method similar to the one used in reference 5 to - 
indicate divergence speed. Often, however, the variation of the rate of 
roll through the speed range is of mare interest than the reversal speed 
alone; if this variation is calculated by use of table I, the reversal 
speed is automatically given, if reversal occurs. Because of the linear 
. pb /27 ^ 3 

dependence of 0-r/B on 1 /ma and q, the values of found 

r (pb/2V) R 


requires some discussion. When F = 1 the values of pc^^F 
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from table I far specified values of stiffness anfi altitude c an "be con- 
verted easily to give results for otter values of stiffness anfl altitude. 

Often a solution is desired for the torsional stiffness required to 
maintain a certain percentage of the rate of roll that would he obtained 
with a rigid wing. This solution is easily effected by a slight rearrange- 
ment of the last several columns of table I , as shown in the table; 

As an example of the type of data obtained from the use of the 
figures and the computation form, figure 12 shows the results of some 
computations for a typical free— flight research model. 


CONCLUDING REMARKS 


An analysis has been presented of the problem of the loss in rolling 
effectiveness at supersonic flight speeds for the case of the rectangular 
wing. The effects of the loss in lift in the various Mach cones were 
considered in the analysis, and equations were obtained for two assumed 
variations of the spanwise distribution of torsional stiffness. 

A computation form and figures were prepared; with the aid of these, 
calculations can be made without the necessity of referring to the details 
of the analysis. 


Langley Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Langley Field, Va., September 21, 19^8 
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APPENDIX 

BOLLING MOMENT BESOTTING EROM PARABOLIC TWIST 


If the pressure distribution over the twisted wing can "be found, 
this pressure distribution may be integrated to give the rolling moment. 
Perhaps the most direct way to find the pressure distribution for the 
particular case in the present paper is -to consider the twisted wing to 
be represented by a distribution of sources and sinks, so that the 
pressure distribution can be found by use of the methods of references 6 
and 7- 

In general, only one wing panel' need be considered (see fig. 13). 

The distribution of vertical velocity along the span of this panel is 

Wg is a constant; on the other panel the distribution 

point located between the Mach cones from the tip and 

(see fig. 13 (a)), the velocity potential on one surface 


0l(x,y) 




pSl 


ft 7 !- 


w, 


_2 

it 


^(^T)) <U dT l - “f 


' UT] 1 ^0 


3 

liy (Jo 




iff(S,Tl) « dT) (Al) 


just WgT] 2 , where 

O 

is — w tj . At any 

center of the wing 
is given by 


where 


f(S,T)) 


(x -t) 2 - 3 2 (y 



- 1/2 


At any point located within the tip Mach cone, the velocity potential is 
given by 



1 6 
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fis 


w_ 


0 2 (x,y) = ^(^y) + — 


nl 


Tl 2 f(ljTl) at dTl 


^2 J 


w, 


P''3 


n£- 


2_ 

ji 


ti as aij 


Js UQ 


(A2) 


where the last two terms represent the effect of the finite tip. At any 
point within the Mach cone from the center of the wing, the velocity 
potential is 


0 3 (**y) = ftjfoy) + 



t] f(l,T]) at at| 


(A3) 


■tfhere the last term represents the effect of the opposite wing panel. 
If a pressure coefficient is defined as 


P 

1 


(Aft) 


where P is the lifting pressure, then 


C 


P 


ft 

V 


(A5) 
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When equations (Al), (A2), and (A3) are substituted in equation (A5) 
and the indicated operations are performed, the pressure distribution 
over the region between the Mach cones from the tip and center of the 
wing becomes 


, _ v 2 2x 2 + 

'Pi Y p3 


(AS) 


the pressure distribution over the part of the wing within the tip Mach 
cone is given by 



JtV 


oob -1 2Py 

P 3 


20s + x 


40y + h0s — 2x 



\[x 2 - p 2 ^2y - 


2 s 



(A7) 


and in the center Mach cone the pressure distribution is 



^2 12y 

nY p 2 


Jx 2 -? 2 ^ + §ej£. + * 2 Bin " 1 

p3 


(AB) 


Wow the rolling moment of one wing panel is simply 


h 

1 


Cpy dy dx 


(A9) 
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where the integration is taken over the area of the panel. The limits 
of integration are shown in figure 13(h). Equation (A9) becomes 


nc 


c/3 


lio 


C_ y dy to + 

3 


n c 


-I 




c. y dy + 


Jo Jx/p 




Cr, y dy to (ALO) 
y 2 


Jo 

P 


These integrations, which, although tedious, can he carried out with the 
help of a table of integrals, give 


h 3\ 

2q_ w 2 


Ip^c - |ps3 c 2 


+ -5_s 2 C 3 _ X Bci 

2k 32 p 


21 c5 
1280 ^2 


The total rolling-^noment coefficient is found by multiplying by 2 (for 
two panels) and dividing by Sb = hs 2 c. After a little simplification, 
the coefficient becomes 


W 2b 2 l 


P°z - T IT 


1 — — + 2. — - — 

PA 3 (pA ) a 


1 1 __ 1 1 
2 (PA) 3 B (PA) A 


(All) 


which is equation (19). 
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TABLE H 

B MD pq AS FUNCTIONS OF M 



Pq far various altitudes 


10,000 ft 20,000 ft 30,000 ft 40,000 ft 


152 

262 

386 

527 

69O 

877 

1,090 

1,331 

1,600 

1,900 

3,925 

6,982 

11,273 


0.458 

821 

564 

377 

244 

.663 

l,4l4 

973 

650 

420 

.831 

2,080 

1,430 

956 

617 

.980 

2,845 

1,956 

1,307 

845 

1.118 

3,726 

2,562 

1,712 

1,106 

1.249 

^,736 

3,256 

2,176 

1,405 

1.375 

5,886 

4,047 

2,705 

1,746 

1.497 

7,184 

4,940 

3,301 

2,132 

I.616 

8,641 

5,942 

3,971 

2,565 

1.732 

10,262 

7,056 

^,715 

3,045 

ro 

8 

0 

21,201 

14,578 

9,739 

6,293 

2.828 

37,700 

25,924 

17,322 

11,188 

3-355 

60,876 

4l,86o 

27,967 

18,067 


50,000 ft 


94 























(a) P = l. 

Figure 2.- Comparison of exact and approximate ' 

b R 

A = 8,- -f = 0 
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Figure 3-- Comparison of "exact” and approximate twist distribution. xhe "emot" results are 
based upon the appro xima te torque distribution shown in figure 2. 
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Altitude , ft 


Figure b.- Variation of dynamic-pressure ratio with altitude at constant 

Mach, number. 




















Mi 








■■■■■■ ■ —■■■■■■■■■■■■■■■a 















